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Two components in charged particle production in heavy-ion collisions. 
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Transverse momentum spectra of charged particle production in heavy-ion collisions are considered 
in terms of a recently introduced Two Component parameterization combining exponential (“soft”) 
and power-law (“hard”) functional forms. The charged hadron densities calculated separately for 
them are plotted versus number of participating nucleons, N par t ■ The obtained dependences are dis¬ 
cussed and the possible link between the two component parameterization introduced by the authors 
and the two component model historically used for the case of heavy-ion collisions is established. 
Next, the variations of the parameters of the introduced approach with the center of mass energy 
and centrality are studied using the available data from RHIC and LHC experiments. The spectra 
shapes are found to show universal dependences on N var t for all investigated collision energies. 


I. INTRODUCTION 


tributions: 


Two-component models have been used in heavy-ion 
phenomenology for a long time. The reason for that is 
that there is no single theoretical approach that can si¬ 
multaneously describe both low-p^ and high-p^ hadron 
production. The main object of study of such models PQ 
is charged particle density, dN c h/dr], which is expected 
to scale with number of participating nucleons, N part , 
or number of binary parton-parton collisions, N co u, for 
“soft” and “hard” regimes of particle production, respec¬ 
tively. Such scaling becomes a subject of various phe¬ 
nomenological discussions - linear scaling with N part is 
expected for “soft” processes, while scaling with N co u 
is expected for the “hard” regime of particle produc¬ 
tion [2Q 

dN ch /drj = A ■ N part + B ■ N coil . (1) 

Recently, another two component approach account¬ 
ing for another aspect of charged particle production - 
transverse momentum spectra d 2 a / dp\dp - has been in¬ 
troduced by the authors [3]. Remarkably, it was also sug¬ 
gested to consider two sources of hadroproduction related 
to “soft” and “hard” regimes, respectively, and therefore 
parameterize transverse momentum spectra by a sum of 
an exponential (Boltzmann-like) and a power-law pp dis¬ 
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1 Also note that N co u oc 
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where Epkin = \/Pt + M 2 ~ M with M equal to the 
produced hadron mass and A e , A,T e ,T, N are the free 
parameters to be determined by fit to the data. 

Moreover, this approach was shown to effectively de¬ 
scribe heavy-ion collision data (?] when the exponen¬ 
tial term of ([2]) is substituted with the well-known Blast- 
Wave formula [5J: 

Pt sinh p \ / niT cosh p \ 
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( 3 ) 

taking into account hydrodynamical expansion of the col¬ 
liding system. In this approach the expanding under the 
pressure in the longitudinal direction system generates 
the transverse flow. The particle distribution is consid¬ 
ered to be Boltzmann again but in the local fluid rest 
frame. In <[3|) p = tanh~ * 1 /3 r and (3 r (r) = /3 S (-^), with f3 s 
standing for the surface velocity, tot = \Jm 2 + p\, Iq 
and I\\ are the modified Bessel functions. 

In [3j it was also shown that an additional power-law 
term is needed to describe the charged hadron spectra 
in central PbPb collisions in the full range of transverse 
momenta. Thus, the experimental data are fitted to the 
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A typical charged particle spectrum as a function of 
transverse momentum measured by the ALICE collabo¬ 
ration in PbPb collisions is shown in figure |T] together 
with the introduced new three component fit function. 


discussed two component approaches ([T]) and 0 can be 
established giving the further insight into the interplay of 
“soft” and “hard” regimes of hadroproduction in heavy- 
ion collisions. 

Let us now have a look how charged hadron densi¬ 
ties dN c h / di] vary with the centrality of heavy-ion colli¬ 
sions. Therefore, to study this scaling we suggest to plot 
charged hadron densities 1 / N evt dN/ dr] / Npart calculated 
separately for “soft” and “hard” parts of 0 as functions 
of N part as shown in the figures [2] and [3j 



FIG. 1. (Color online) Central lead-lead collisions $3]: the 
red (dashed) line shows the hydrodynamic term 0 and the 
green (solid) and blue (dash-dot) lines - two power-law terms 
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While the high —pt tail (shown with the green line 
in figure |T|) has the same slope parameter N as in 
pp— collisions [40 and can be rather well reproduced 
in Monte Carlo generators, the low-pr (“soft” part) and 
mid-pr (“hard”) parts (shown with red and blue lines in 
figure 0) of the spectra have only phenomenological de¬ 
scriptions. Therefore, in this paper we suggest to look 
how these two contributions and their parameter values 
vary with centrality and center of mass energy to better 
understand the dynamics of charged hadron production 
in heavy-ion collisions. 


II. CHARGED HADRON DENSITIES 

First of all, recall that integrating the transverse mo¬ 
mentum spectrum d 2 a/dp^drj by dp^ from 0 to the up¬ 
per kinematic limit, one gets the charged particle den¬ 
sities dcr/<i ?7 0 Therefore, the possible link between the 



FIG. 2. (Color online) The charged particle densities divided 
over the number of participating nucleons N par t obtained for 
the mid-pr part (blue line in the figure 0 of the spectra as a 
function of N part . The lines show the scaling of this quantity 
as N p the same for RHIC [7] and LHC [6] data. 


As one can notice, in figure [2] charged hadron densi¬ 
ties obtained for the power-law term contribution show 
an universal scaling with Np art (a = 1.14) independent 
on the collision energy, with a > 1 indicating that the 
mid-pr power-law terms of §■ @ are related to some 
“hard” regime as postulated in our model. The value of 
a < 1.33 which might be expected for the purely “hard” 
regime of hadroproductiorjj is explained by the quench¬ 
ing of mini-jets in heavy-ion collisions. The rescattering 
of the original jet diminishes the mean pt of the jet but 
enlarges the number of the final mini-jets, that leads to a 
larger particle density dN/dr]. This is the origin of an ad¬ 
ditional (in comparison with 0) power-law term in (4) 
which contribution is shown by the blue line in figure 1 
At first sight the exponential part caused in pp- case 
just by the incoming valence quarks [S] should be pro¬ 
portional to Np ar t ■ This is more or less consistent with 
the behavior observed in figure [3] for the RHIC energies. 
However, at a higher, LHC energy, the energy density of 
the partons increases and the number of final state rescat¬ 
terings becomes so large that the secondaries produced 
in these rescatterings start to thermalize. In our fit this 
contribution is described by the first term of 0. By this 


2 This observation is also supported by the fact that the nuclear 
modification factor Raa comes to a constant value for the high- 
Pt region. 

3 The rjiab = 0 region is considered. 
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FIG. 3. (Color online) The charged particle densities divided 
over the number of participating nucleons N par t obtained for 
the “thermal” part (red line in the figure [l]) of the spectra as 
a function of N par t- 


reason the ’thermal’ part of the particle density dN/dr] 
increases with N part , especially at the LHC energy. 

To perform the quantitative test we suggest to calcu¬ 
late the total transverse energy EdEx/dp for the “ther¬ 
mal” part of the spectra. This quantity divided over the 
number of participants is shown in figure [4j One can see 
that the total transverse energy EdEx/dp has the same 
scaling with N p ^}j. for both RHIC and LHC data which 
confirms the picture described above. 

Remarkably, the observed scaling EdEx/dr] oc ss 

N/ a /i correlates with the fact that naively EdEx/dp 
should be proportional to the total energy stored in the 
colliding system. This observation further supports the 
idea that the two component model indeed reveal the 
underlying hadroproduction dynamics in heavy-ion colli¬ 
sions and not just a useful ansatz as was argued in [9]. 



FIG. 4. (Color online) The total transverse energy for “soft” 
particle production (red line in the figure [I]) divided over the 
number of participating nucleons N part as a function of N part . 
The lines show the scaling of this quantity as JV p art the same 
for RHIC and LHC data. 


III. PARAMETERS 

Let us now have a look how the values of the param¬ 
eters T and N of the power-law term of Q (blue line 
in figure |T|) depend on the energy and the centrality of 
a collisions. In [3j it was suggested that the mini-jets 
travelling through the dense medium should loose their 
energy because of the multiple rescatterings. This ef¬ 
fect, in turn, will result in higher values of T and N 
parameters (of the mid-pr term), as expected for higher 
thermalization. Here we consider pp— collisions as a ref¬ 
erence point for the case of heavy-ion collisions to study 
the modification of the spectra with centrality. On the 
other hand, the most universal parameter characteriz¬ 
ing the centrality of a heavy-ion collision is N part or the 
number of participating nucleons related to the impact 
parameter. Therefore, we propose to look at the varia¬ 
tion of the ratios T/T pp and N/N pp obtained from the fit 
to the experimental data as a function of N part , where 
T pp and N pp stand for the parameter values obtained 
for pp—collisions at the same c.rn.s. energy per nucleon 
Vs/N. 



FIG. 5. (Color online) The dependence of the ratio of the 
parameters N and N pp of the two-component model on N par t- 
The N is obtained for heavy-ion collisions and N pp for pp- 
collisions at the same c.m.s. energy. The red (solid) line 
shows the dependence 

As shown in figures [5]|6] indeed, both of the parame¬ 
ters, T and N , grow with N part as it was natural from our 
naive expectations. Using the data from the ALICE [6], 
ATLAS pH] and PHENIX [7j experiments their behavior 
can be parameterized in the following way as a function 
of N part : 


N/N pp = 0.02 • N° p ™ + 0.97, (5) 

T/T pp = 0.07 • N pap ® + 0.91. (6) 
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FIG. 6. (Color online) The dependence of the ratio of the 
parameters T and T vv of the two-component model on N par t ■ 
The T is obtained for heavy-ion collisions and T pp for pp- 
collisions at the same c.m.s. energy. The red (solid) line 
shows the dependence |[6j. 


Remarkably, in figures [5]|6] the points corresponding to 
pp-collisions are nicely placed on these trends ex¬ 

tracted from the heavy-ion data. Moreover, from these 
figures [5]|6] one can notice a universal scaling of the pa¬ 
rameters T and N with N part to be the same both for 
RHIC and LHC data. This supports the suggested pic¬ 
ture of hadroproduction in heavy-ion collisions and there¬ 
fore one can conclude that the change of the spectra 
shape with centrality is indeed related to the higher sup¬ 


pression of the mini-jets produced in a collision. 


IV. CONCLUSION 

In conclusion, transverse momentum spectra in heavy 
ion collisions have been considered using the two compo¬ 
nent parameterization. Scaling of exponential and power- 
law contributions to the charged hadron spectra with 
centrality in heavy-ion collisions were studied separately. 
The charged hadron densities for “hard” part turned out 
to have the same scaling for RHIC and LHC data with 
a > 1 indicating the “hard” regime of hadroproduction. 
For the “soft” part the deviation from awl was found, 
especially for the LHC, and is explained by large initial 
parton densities which strongly enhance the rescattering 
of the mini-jets leading to the thermalization of the part 
of secondaries, that is to the additional contribution to 
the “soft” term. Finally, the variations of the parame¬ 
ters obtained from the fit have been studied as a function 
of Npart and c.m.s. energy. The same universal depen¬ 
dence of the parameters of the model, T and N, for RHIC 
and LHC data was found giving further insight into the 
hadroproduction dynamics in heavy-ion collisions. 
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